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The coefficient of the logarithmic term in the entropy on even spheres 
is re-computed by the local technique of integrating the finite tem- 
perature energy density up to the horizon on static <i-dimensional de 
Sitter space and thence finding the entropy by thermodynamics. Nu- 
meric evaluation yields the known answer i.e. (minus) the conformal 
anomaly on the (i-sphere. The de Sitter quantities are obtained by 
conformal transformation of the Rindler ones, themselves obtained, for 
convenience, from those around a cosmic string. The expressions are 
given in terms of generalised Bernoulli polynomials for which an iden- 
tity is derived. The arising spherical conformal anomaly is discussed 
and a formula is given for it for Branson's higher GJMS Laplacian, 
P 2 k, as an oscillating polynomial in the level, k. 
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1. Introduction. 

In a previous work, [1], I have computed the universal logarithmic term that 
occurs in the expansion of the entanglement entropy associated with even spheres 
in the special case when the conically deformed Euclidean space-time used in the 
entropy recipe is the cyclically factored sphere, S d /Z q . Then the spatial surface is 
a (d — 2)-sphere (with vanishing extrinsic curvature). 

In this note I return to the same problem but the method this time uses the 
thermodynamic way of introducing a conical singularity, [2], and harks back, in de- 
tail, to a computation of the geometric entropy on de Sitter space, [3] , by straight- 
forward integration of the energy density off shell and then picking out the log term 
in the entropy before going back on shell. 

This off shell procedure is well known and is used mostly in connection with 
the black hole geometry. Fursaev and Miele, [4] , discuss the specific application to 
de Sitter space, relevant to the present situation. See also De Nardo et al, [5]. 

A recent paper by Marolf et al also sets up thermodynamics on static de Sitter 
by integrating the energy density holographically generated. 

Calculations of spherical entanglement entropy have been performed by Casini 
and Huerta, [6], Solodukhin [7,8] and Lohmeyer et al., [9]. 

2. The log coefficient in four dimensions 

The geometry I treat is de Sitter space and the method, and formulae, are 
to be found in [3]. The novelty here is that the expansion away from the horizon, 
which there stopped at the most divergent term, is here taken far enough to give the 
logarithm. This is somewhat elementary but the different approach makes contact 
with the more geometric and holographic methods. Similar, but not identical, 
calculations have been done over the past fifteen years. 

I start by writing the de Sitter metric in the static form, 

ds 2 = -^^y {Z 2 d{t/af - dZ 2 ) - a 2 (j^) (dO 2 + sin 2 # 2 ) , (1) 

where Z is related to the usual static radius r by Z 2 = (a — r)/(a + r). For conve- 
nience I set a, the radius of the de Sitter sphere, to unity. There is a slight com- 
putational advantage in using Z rather than r, and the metric exhibits a Rindler- 
like part. The part of the manifold I am interested in corresponds to the range, 
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< Z < 1, the lower limit being the location of the horizon where the field the- 
ory averages diverge. In this section I use a lower cut-off, Z > e/2, to define the 
global quantities. The basic idea is to compute the energy density, (Tq)^, at a finite 
temperature, T = 1/(3. Then integrate this over the region indicated to give the 
internal energy, E(f3), and next find the entropy from, 

S{0) = 00(0) - f E'(p) dp = J% dE\p) , (2) 

where E'(/3) is the finite temperature correction to the energy. 

The last act is to go on shell by setting /3 = 2tt, the Gibbons-Hawking temper- 
ature associated with the horizon. At this temperature, static de Sitter turns into 
zero temperature de Sitter, [10], the state being the de Sitter invariant 'Euclidean' 
vacuum because the corresponding Green functions are related by the thermalisa- 
tion, [11] 2 . See also Laflamme, [12]. 

Expressions for (Tq)^ in four dimensions were given in [3]. There are various 
means of finding them. They can be calculated directly in de Sitter space from the 
re-periodised Green function, [11]. They can be found by conformally transforming 
the expressions in Rindler space, which can themselves be obtained from those 
around a cosmic string by transposition of coordinates. This last is the path I 
will follow here because the cosmic string results have already been calculated in d 
dimensions, [13]. I give a few subsidiary details. It will be noticed that B in the 
present situation plays the role that q did in [1] . 

For convenience I repeat the formulae on Rindler space with metric, 

ds 2 = Z 2 d t 2 - dZ 2 - dx 2 - dy 2 , (3) 

(T )^^(5 2 -l)(5 2 + l), j = 

' (£ 2 -l)(7B 2 + 17), j = 1/2 (4) 



1920tt 2 Z 4 



1 (B 2 -1)(B 2 + 11), j = l, 



240tt 2 Z 4 
where B = 2n//3. 

These quantities vanish on shell, i.e. when B = 1, which is a consequence 
of the subtraction of the Minkowski values. In fact the last terms in (4), the zero 
temperature values, do not contribute to the entropy, (2), as is easily seen, and we 
could simply concentrate on the remaining finite temperature corrections. 



2 However see [4]. 
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The transition to de Sitter space uses the Bunch-Davies-Brown-Cassidy (BD- 
BC) relation between the averaged energy momentum tensor densities on confor- 
mally flat manifolds developed in [10], where the original references are given. It 
is not necessary to give this relation in detail. It connects tensor densities under 
conformal transformations in the expected manner, but there is an additional piece 
resulting from the integration of the conformal anomaly. For spaces conformal to 
Rindler, as here, with static conformal factors, the relation also holds for finite 
temperature quantities. The addition, however, is independent of temperature and 
can be ignored when calculating the entropy. The finite temperature corrections 
are thus homogeneously related by 3 , 

V9(T^ = V9 R (Tn, (5) 

where R refers to Rindler (the expressions (4)) and the dash signifies the finite 
temperature corrections. 

The relation of the coordinates implied in (5) is explained in [10]. The upshot 
is that the finite temperature corrections in static de Sitter can be obtained from 
(4) by the replacement, 

■1 ^ 1 + ^ 2 ) 4 (6) 
The total de Sitter thermal energy is the integral, 

E' J (P)= f V9 dZ d9 dcf> <T % 

Je/2 

{p d (B) - Pjjoo)) r 1 (i-z 2 ) 2 

167T J £/2 Z* 



dZ, 



where Pj(B) is the polynomial in (4), j indicating the spin i.e. 



{T R, 0) = 



16tt 2 Z 4 ' 



The constant, zero de Sitter temperature, term, P-(oo), will not contribute to the 
entropy when (2) is constructed and hence can be dropped. 
It then is easy to show that, 



Sj(B)= I -dPj(B) j e/2 ^T^dZ, (7) 



This is already evident for global quantities from the considerations in [14] 
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which, on shell, gives the simple expression, (using -P, (l) = 0), 



p.(B) f 1 (1-Z 2 ) 2 



I concentrate on the loge term as the others are not universal. Its coefficient 
is the number, 

1 f 1 dB 



which is easily computed. I find 



1 r 11 r - 16 

90' Cl / 2 ~180' Gl "45 



and I should comment on these values. 

The first two agree with (minus) the accepted values of the conformal anomaly 
on the four sphere if I note that I am using a Weyl neutrino, as in [10]. However 
the spin one value differs from the standard Maxwell conformal anomaly, —62/90, 
e.g. [15]. In [16], a relevant heat-kernel coefficient (that for transverse vectors) is 
given as —2/45. 



3. Extension to higher dimensions 

I now extend the procedure to d dimensions, but only for conformal scalars. The 
two factors in (8) extend separately. This is similar to the analysis of Solodukhin, 
[8] , and Casini and Huerta, [6] , where the formulae, following the geometry, likewise 
split into two. 

The second integral in (8) is geometric and the extension to d dimensions 
consists of replacing the two-sphere metric in (1) by that for a (d — 2) -sphere. This 
trivially alters the \/Jg~\ measure. 

The first integral is field-theoretic. I will again obtain it from the expressions 
for (Tq) around a cosmic string in d dimensions, [17], via the Rindler values which 
are related by, 

(?o) Rindler = — — 1) C^o) string ■> 

using tracelessness of the local (T"). 
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The general structure is now, 

o _ P(d,B) 

\ /Rindler 2 d TX d / 2 Z d ' 

The BD-BC relation, (5), is still valid for the finite temperature corrections, 
which is all I need. It is fortunate the conformal anomaly contribution is not required 
since its <i-dimensional form is unknown. 

After the Rindler to de Sitter replacement, cf (6), 

j_ (i + zy 

Z d ^ 2 d Z d ' 1 ' 

the expression for the entropy on shell becomes, as before, 

\ sd ~ 2 \ /^^S) f 1 (l-Z 2 ) d ~ 2 



SAP) 



Picking out the log term in (10) I find its coefficient to be, 
C{d) ~ 2^-3^/2-1 [d/2 -l)J dB 



2 d - 2 T(d/2) J B 



(11) 



Perhaps a more fundamental quantity is the free energy and I just write down 
its off shell finite temperature correction, 

F ^ B ) ~ 2 2 d - H a/2-2 J dB ^2— J £/2 zd -i dZ - ( 12 ) 

The polynomial, P(d, B), is actually the combination of two polynomials (see 
below). However, for the moment, I just present some explicit examples which can 
be read off from the formulae in [13], especially equns. (10) and (ll). 4 

I list a few of the P(d, B), 

P(6, B) = ^—{B 2 - 1)(10S 4 + 31S 2 + 31) 
v ' 945 

P(8,B) = ^^(S 2 -1)(7S 2 + 17)(3S 4 + 10S 2 + 17) (13) 
P(10, B) = -^^(B 2 ~ 1)(30S 8 + 261B 6 + 1031S 4 + 2219B 2 + 2219) . 



The expressions in [13] are for complex scalars. Hence there is a factor of two difference. 
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The polynomials do not contain a B 2 term, which is generally true and a 
consequence of conformal invariance (for scalars). 

Compact expressions can be given in terms of generalised Bernoulli polyno- 
mials, as briefly mentioned in [13]. Algebra produces the explicit formula in d 
dimensions as the combination of two polynomials, 

m B) = _ ?qzWZ?) U t » mB , i) + «£z*L*& m - mb, d) , 

(14) 

and it is this combination that accounts for the cancellation of the B 2 terms. Each 
term in (14) vanishes at B = 1 by virtue of a standard property of Bernoulli 
polynomials, [18], [19]. 

Substitution into (11) then produces, on shell, the log coefficient, 

C < d > = f § (^'Wm 1) + W 2 - l\B, 1)) , 

(15) 

whose numerical evaluation yields (minus) the scalar conformal anomalies on the en- 
sphere, as do the calculations mentioned earlier. The bare evaluation itself gives no 
clue as to why this should be so. In our previous approach, [1], the analytical reason 
was that the derivative of the conformal anomaly on the orbifolded <i-sphere was 
zero for the non-singular sphere, i.e. on shell, and I make some relevant remarks 
in the next section. 

The conformal anomaly values are available in the old literature and, for con- 
venience, and interest, I give a few remarks and specifics in the next section. 

4. Conformal anomaly on spheres 

As has been used above, the conformal anomaly is the constant term in the 
heat-kernel expansion. 5 There are a number of ways of finding this, based on the 
standard spectral information for the (conformal) Laplacian on <i-spheres (which 
can be traced back to George Green in the 1830s). 

There is no need to elaborate, as the work of Copeland and Toms, [20], cov- 
ers the necessary ground. One technique uses the values of the corresponding £- 
function, the required analytic continuation of which was done by expressing it as 
a series of Hurwitz ^-functions, which have known properties. Birmingham, [21], 

5 The trace anomaly is a constant times the conformal anomaly. 
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employs a similar approach. A list of the numerical results can be found in [20] for 
scalars and spinors. 

Among other works, I mention the interesting one by Cappelli and D'Appolonio, 
[22], which employs basically the same method, only slightly more streamlined. It 
also lists the values. 

A different organisation of the spectral information was given in [23,24] and 
resulted in the compact expression for the heat-kernel coefficients (the conformal 
anomaly) on a (i-sphere, 

C d/ 2 = ^ (Bf\d/2 - 1|1) + Bf\d /2\1)^ , (16) 

which rapidly evaluates to the known numbers. Spelling out some calculational 
details, this evaluation can be done using the integral relation, [18], p. 171, 



l n - 1 \x\l)= f dtBi n \x + t\l), 
Jo 



B 

which produces, from (16), 

C d/2 = I jf 1 dt [Bf +1 \d/2 + 1 11) + Bf +1 \d/2 - 1 11)) , (17) 

with the product form, {e.g. [18] p. 186, [25], §8), 

B { f +1 \x\ 1) = (x - l)(x - 2) . . . (x - d) . 
This yields an expression identical to that given by Diaz, [26], 

d/2-l 



c*/2 = *^jr [ dt II (* 2 -* 2 )' 



(18) 

i=0 

which he obtains on the basis of a holographic technique as the simplest example 
of a more general construction of the conformal anomaly associated with Branson's 
GJMS higher Laplacian operators, P 2fc , on the sphere, e.g. Branson, [27], Graham, 
[28], Gover [29], Juhl, [30]. 

It is possible to show 6 , using the approach of [24], that the general result at 
level k is, 

C d/2 (k) = ^f£ (B d d \d/2 + j |1) + B d d \d/2+j + 1 11)) 

(19) 

2(-i)«*/2 /■* d/ £r\. 2 2 , 



6 The details will be presented at another time. 
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which evaluates to a polynomial in k, for a given d; for example, 



C 2 (k) = ^k 2 (3k 2 -5) 

C 3 (k) = k 2 (3k 4 - 21k 2 + 28) 

C M = tAt fc2 ( 5A;6 " 90A;4 + 441fc2 " 54 °) 

5.9! {20) 

C 5 (k) = k 2 (3k 8 - HOfc 6 + 1287A; 4 - 5412A; 2 + 6336) 

C 6 (fe) = /c 2 (105/c 10 - 6825/c 8 + 155155/c 6 - 1490775/c 4 

105. 13! 

+ 5753748/c 2 - 6552000) . 

Evaluation gives agreement with the numbers found by Diaz using holography, [26] , 
eqn.(7.1) and Table 3. 

The extrema of these polynomials lie in the region \k\ < d/2 and are located, 
for positive k, closely below the integers which fact corresponds to the alternating 
signs at the 'physical' values, k = 1, . . . , d/2. 

In this approach, the operator existence condition on the level, k < d/2 (for d 
even), is equivalent to the absence of a zero mode in the spectrum of P 2 fc- (Actually 
for the Neumann modes on the hemisphere.) It may be possible to relax this 
condition. 

It is worth noting that (17), (18) and (19) do not involve the Bernoulli numbers, 
unlike some previously mentioned computations. 

Returning to the simplest case, k = 1, equating (16) and (15) gives an identity 
between generalised Bernoulli polynomials. In fact there is the off shell identity 



i (Bf\d/2 - 1 1 B, 1) + Bf\d/2 | B, 1)) = 

W -Vf§ W2|5, 1) + d -^B^\d/2 - 1)) . 

(21) 

The left-hand side is d\ times the conformal anomaly on the orbifolded en- 
sphere (or periodic lune), S d /Z B , discussed in [1] where it was proved that it had 
an extremum when B = 1, i.e. on shell. This can be confirmed here immediately 
since (each term in) the integrand vanishes on shell. Reference to (12) shows that 
this corresponds to an extremum for (3 F, in thermodynamic language. 
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By applying Barnes' derivative formula, 

+ B\B,B,1), 

(22) 



B, 1) = -± B[ d+1 \a + B\ 5, 5, 1) , 



= "(-I)' ^ (d - 1 - a + B\ B, B, 1) , 

to (21), I obtain the off shell identity (d is even), 

B^\d/2 | 5, l) + ^_L^ 5^_- 1 )(rf/2 - 1 1 5, 1) 



1 | (23) 



which is not obvious without further manipulation. It furnishes a simpler formula 
for the scalar conformal polynomials, P{d,B), (14), 

P(d, B) = ~ 2d " 1 ^ ( 1 j! 2 i r ) f /2) Bf + »(d/2 + B\ 5, 5, 1) . (24) 



5. Comments 

The calculation of the entropy in section 4 has similarities with others in this 
topic, e.g. [8,6]. However, in these, the entropy is usually computed from the global 
effective action while I have employed a local method which involves integrating 
the thermal average of the canonical energy obtained, e.g., from the coincidence 
limit of a bilinear operator acting on the thermal Green function. This gives the 
total internal energy from which the free energy can be determined and thence the 
entropy. 

The numerical evaluation is algorithmic and offers no obvious reason why the 
log coefficient is (minus) the conformal anomaly, although, being conformally in- 
variant, there is not much else it could be. In (15) one has a rather complicated 
formula for calculating the conformal anomaly on spheres. 

I have taken a simple minded approach to the conical singularity and have inte- 
grated the local energy density using a basic cut-off, resorting neither to smoothing 
the geometry nor to integrating against a test function. These procedures might 
appear more rigorous, but have their own ambiguities. I believe the present method 
is sufficient to achieve its aim. 
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I also note that the formulae encountered are closely related to expressions 
occurring in the holographic calculations of Ryu and Takayanagi of the entangle- 
ment entropy for a disk, [31] §7.1.2, and also to those of Graham, [32], on volume 
renormalisation. 

The finite temperature averaged stress tensor on static de Sitter that I have 
used does not agree with that derived by Marolf et al, [33], on the basis of a dS/CFT 
duality. 
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